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Abstract. We present an 0(n log n) time algorithm for the (weighted) 
fc-center problem of n points on a real line. We show that the problem 
has an J?(nlogn) time lower bound, and thus our algorithm is optimal. 
We also show that the problem is solvable in 0(n) time in certain spe- 
cial cases. Our techniques involve developing efficient data structures for 
processing queries that find a lowest point in the common intersection of 
a certain subset of half-planes. This subproblem is interesting in its own 
right and our solution for it may find other applications. 



1 Introduction 

We study the weighted fc-center problem for a set of n points on a real line. 
Let P — {pi, p2, ■ ■ ■ ,p n } be a set of n points on a real line L. For each i 
with 1 < i < n, the point pi £ P has a weight w(pi) > 0. For a point p on 
L, denote by L(p) the coordinate of p on L, which we also refer to as the In- 
coordinate of p. For two points p and q on L, let d(p,q) = \L(p) — L(q)\ be 
the distance between p and q. Further, for a set F = {fx, / 2 , . . . , fk} of points 
and a point q on L, define d{q,F) = d(F,q) = mini<y<fc d(q, fj). Given P 
and an integer k > 0, the weighted one- dimensional k-center problem seeks 
to determine a set F = {fi, f2, ■ ■ ■ , fk} of k points on L such that the value 
tp(P,F) — max Pi £p(w(pi) ■ d{pi,F)) is minimized. We use IDkC 'enter to denote 
this problem. Also, the points in F are called centers, and the points in P are 
called demand points. 

The unweighted version of lDkCenter is the case where all points have the 
same weight. If F C P is required, then the case is called the discrete version. 
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Although many variants of the /c-center problem are NP-hard |H3ll2ll5ll7j . 
some special cases are solvable in polynomial time. Megiddo and Tamir [T7] pre- 
sented an 0(n log 2 n log log n) time algorithm for the unweighted fc-center prob- 
lem on a tree of n nodes, and the running time can be reduced to 0(n log 2 n) 
by applying Cole's parametric search [7]. Later, Frederickson [TT| gave a linear 
time algorithm for the unweighted fc-center problem on a tree. Jeger and Kariv 
|14j gave an O(fcnlogn) time algorithm for the weighted fc-center problem on a 
tree. For the weighted fc-center problem on a real line (i.e., the problem Wk- 
C enter), Bhattacharya and Shi [5] recently proposed an algorithm with a time 
bound linear in n but exponential in k. In addition, the discrete weighted k- 
center problem on a tree is solvable in 0(n log 2 n) time |18j and the discrete 
unweighted fc-center problem on a tree is solvable in 0(n) time Tamir [19] 
indicated that the discrete weighted WkCenter can be solved in 0(n\ogn) time 
by using the procedure given in the last section of [T8]. Below, unless otherwise 
stated, WkCenter always refers to the weighted version. 

In this paper, we solve the lDkCenter problem in O(nlogn) time, which is 
faster than the 0(kn\ogn) time result [T3] when they are applied to WkCenter. 
Our result is also better than the solution in [5] for large k since the running 
time in [5] is exponential in k. In fact, we can show that the problem has an 
Q{n\ogn) time lower bound, and therefore, our 0{n\ogn) time algorithm is 
optimal. Further, if all points in P are given sorted on L and their weights are 
also sorted, we solve WkCenter in 0(n + k 2 log 2 ? log n log log n) time, which is 
in favor of small k. For example, if k — 0(n 1 / 2 ~ e ) for any e > (which is true 
in many applications), our algorithm runs in O(n) time. 

In addition, our techniques also yield an efficient data structure for processing 
queries for finding a lowest point in the common intersection of a certain subset 
of half-planes, which we call the 2-D sublist LP queries. Since the 2-D sublist 
LP query is a basic geometric problem, our data structure may be interesting in 
its own right. 



1.1 An Overview of Our Approach 

We model the WkCenter problem as a problem of approximating a set of 
weighted points by a step function in the plane [5] . Two algorithms were given in 
[5] for this point approximation problem with the time bounds of 0(n log 2 n) and 
0(n log n + k 2 log 2 ? log 2 n), respectively. Consequently, the WkCenter problem 
can be solved in 0(min{n log 2 n, n log n + k 2 log 2 ? log 2 n}) time. 

However, the WkCenter problem has some special properties that allow us 
to develop faster solutions. Specifically, after the geometric transformations, a 
key component to solving the problem is the following 2-D sublist LP query 
problem: Given a set of n upper half-planes, H = {/ii, /12, . . . , h n }, in the plane, 
for each query q(i,j) (1 < i < j < n), compute a lowest point p* in the common 
intersection of all half-planes in Hij = {h t \ i < t < j}. A data structure 
was proposed in [5] for this problem, which can be built in O(nlogn) time and 
answers each query in 0(log 2 n) time. On the WkCenter problem, we observe 
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that the input half-plane set H has a special property that the intersections 
between the x-axis and the bounding lines of the half-planes are ordered from left 
to right according to the half-plane indices in H. Exploiting this special property 
and using the compact interval trees |13j . we design a new data structure for 
this special case of the 2-D sublist LP queries, which can be built in O(nlogn) 
time and can answer each query in 0(log n) time. This new data structure allows 
us to solve the IDkCenter problem in O(nlogn) time. Since the 2-D sublist LP 
query problem is a very basic problem, our new data structure may find other 
applications as well. 

We should mention that very recently, by using Cole's parametric search 
[7], Fournier and Vigneron [10] gave an 0(n\ogn) time algorithm for the above 
point approximation problem [5]. Thus, by combining our problem modeling 
and the algorithm in jTU], the IDkCenter problem can be solved in O(nlogn) 
time. However, as pointed out in |10) . the parametric search approach in [7] is 
quite complicated and involves large constants, and thus the algorithm in [10] 
is mainly of theoretical interest. In contrast, our approach is much simpler and 
more practical. 

In the following, we present the high-level scheme of our algorithm in Section 
[5] In Section [31 we model our problem as the 2-D sublist LP queries and present 
our data structure. Section 2] concludes the paper and discusses the lower bound 
of the problem IDkCenter. 

For simplicity of discussion, we make a general position assumption that no 
two points in P are at the same position on L. We also assume the weight of 
each point in P is positive and finite. These assumptions are only for ease of 
exposition and our algorithms can be easily extended to the general case. 

2 The Algorithmic Scheme 

In this section, we discuss the high-level framework of our algorithm. As pointed 
out in [2], it is possible that there is more than one optimal solution for the 
IDkCenter problem. Our algorithm focuses on finding one optimal solution. 

2.1 Preliminaries 

For any two points p and q on L with L(p) < L(q) (recall that L{p) is the 
coordinate of p on L, and similarly for L(q)) 7 denote by [p, q] the (closed) interval 
of L between p and q. 

We first sort all points of P from left to right on L. Without loss of generality 
(WLOG), let {pi,P2, ■ ■ ■ ,Pn} be the sorted order in increasing coordinates on L. 
For any two points Pi,Pj 6 P with i < j, denote by the interval [pi,Pj]. 

Let ip* be the value of ip(P, F) for an optimal solution of IDkCenter. Suppose F 
is the center set in an optimal solution; for a demand point p £ P and a center 
/ S F, if (w(p)-d(f,p)) < then we say that p can be served by /. It is easy to 
see that there is an optimal solution F such that each center of F is in [pi , p n ] . 
Further, as discussed in [2], there is an optimal solution F such that the points 
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Fig. 1. Approximating a set of points by a step function (i.e., the three horizontal line 
segments) . 

of P are partitioned into k intervals 7(1, ii), + l,^), ■ • ■ ,I(ik—i + l>n) by 
integers iq + 1 = 1 < i\ < i 2 < • • • < ik-i < n — ik, each interval + 1, ij) 

contains exactly one center fj S F, and for each point p G P PI F\ij-i + 1, ij), 
(w(p) ■ d(fj,p)) < -0* holds. In other words, each center of F serves a subset of 
consecutive demand points in P. 

For any two integers i and j with 1 < i < j < n, denote by Py the subset of 
points of P in the interval i.e., Py = {pi,pi + \, . . . ,pj} (Py = {p{\ for i = 

j). Consider the following weighted 1-center problem: Find a single center (i.e., a 
point) / in the interval I[i, j] such that the value of ip(Pij > f) = rnax Pt e p 4 . (iy (p t ) • 
d(pt, /)) is minimized. Let a{i,j) denote the minimum value of ip(Pij, f) for this 
weighted 1-center problem. 

For solving the WkCenter problem, our strategy is to determine k — 1 integers 
1 < i-y < i 2 < • ■ • < < n such that the value of max{a(l, ii), a{i\ + 

1, i 2 ), ■■• , a(ifc_i + 1, n)} is minimized and this minimized value is "0*- Note that 
in the above formulation, for each value a(i,j), exact one center is determined 
in the interval I(i,j). To solve this problem, we reduce it to a planar weighted 
point approximation problem [5 in the next subsection. 

2.2 The Reduction to the Planar Weighted Point Approximation 
Problem 

We first review the planar weighted point approximation problem and then show 
our problem reduction. 

Let P' = {p'i,P2, ■ ■ ■ iPn} be a point set in the plane with p'- = (xi,yi), and 
each point p\ be associated with a weight w{p' i ) > 0. Assume the points in P' 
are ordered increasingly by their x-coordinates. Suppose g is a step function 
(i.e., a peicewise constant function, e.g., see Fig.[TJ which we use to approximate 
the points of P 1 (in other words, we fit the step function g to the point set 
P'). The weighted vertical distance between any point p[ £ P' and g is defined 
as d w (p' i ,g) ~ w(p' i ) ■ — g{xi)\ (see Fig. [TJ. The approximation error of g, 
denoted by e(P',g), is denned as max p ' e p/ d w (p' i ,g). The size of g is the number 
of its horizontal line segments. Given an integer k > 0, the point approximation 
problem seeks a step function g to approximate the points of P' such that the 
size of g is at most k and the error e(P,g) is minimized. In [5], this problem is 
referred to as the weighted step function min-e problem, denoted by WSF. Here 
we also use WSF to denote this problem. 
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We now show that the WkCenter problem can be reduced to WSF. Con- 
sequently, WSF algorithms can be used to solve WkCenter. Indeed, consider 
the demand point set P = {p\,p2, ■ ■ ■ ,p n } for IDk Center with the points or- 
dered increasingly by their coordinates on L. For each demand point pi G P, 
1 < i < n, we create a point p\ = (i,L(pi)) in a 2-D Euclidean plane R 2 (i.e, 
the ^-coordinate of p\ in R 2 is the index i and its y-coordinate y\ in R 2 is the 
coordinate of Pi on L), and let the weight of p\ be that of pi (i.e., w(p' i ) = w{pi)). 
Let P' be the set of n weighted points thus created in R 2 . The next lemma states 
the relation between WkCenter and the reduced instance of WSF. 

Lemma 1. An optimal solution OPTp> for WSF on P' in R 2 corresponds to 
an optimal solution OPT p for WkCenter on P. Further, once having OPT pi, 
OPT p can be obtained in 0{n) time. 

Proof: For any two integers i and j with 1 < i < j < n, let P(j = {p' i ,p' i+1 , . . . ,p'j} 
(Py = {p'j} for i — j). Consider the following problem: Find a value Y for one 
single horizontal line segment with Y as its y-coordinate such that the value of 
d w {Plj,Y) = maxpjgjw (wO 5 *) • \y t —Y\) is minimized (where y t = L(p t )). Let 
a'(i,j) denote the minimized value of d w (P!j,Y). 

Let e* be the approximation error of an optimal solution for WSF on P 1 . It 
is easy to see that computing an optimal solution for WSF on P' is equivalent 
to determining k — 1 integers 1 < i\ < %2 < ■ • ■ < ife-i < n such that the 
value of max{a'(l, a'(i\ + 1, 1%), . . . , a'(ik-i + 1,")} is minimized and this 
minimized value is e* . According to the way that we create the point set P' from 
the demand point set P, each value a'(i,j) is exactly equal to the value a(i,j), 
which is the minimized value of ip(Pij,f) for the weighted 1-center problem 
on the demand point subset Py by determining the value of /. Further, we 
have shown that to find an optimal solution for WkCenter on P, it suffices to 
determine k — 1 integers 1 < i± < %2 < • • • < ik—i < n such that the value of 
max{a(l, a(ii + • • ■ ><*(ifc_i + ^-i n )} is minimized and the minimized 

value is tp*. 

The above discussion shows that to find an optimal solution for WkCenter 
on P, it suffices to find an optimal solution for WSF on P'; further, ip* = e*. 
Given an optimal solution OPT pi for WSF on P', below we show how to obtain 
an optimal solution OPT p for WkCenter on P from OPTpi in linear time. 

Note that OPT pi is a step function with k steps (i.e., horizontal line seg- 
ments). Let io = and ik — n. For each 1 < j < k, suppose the j-th step of 
OPT pi has a y-coordinate y- 7 and covers the points of P' from p' ijl+ i to p\., 
i.e., for each point p' t 6 P 1 with +1 < t < ij, the vertical line through p' t 
intersects the j-th horizontal segment of OPT pi . We obtain OPT p for WkCenter 
on P as follows. For each 1 < j < k, the j-th center fj is put at the position yi 
on L (i.e., L(fj) = J/- 7 '), which serves the demand points of P from pu_ x +i to . 

Thus, once OPT pi is available, OPT p can be obtained in O(n) time. □ 

Based on Lemma [U to compute a set F oik centers for P to minimize the 
value ip(P,F), it suffices to solve the corresponding WSF problem on P' and k. 
Specifically, after an optimal step function g for P' is obtained, each horizontal 
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segment of g defines a center on L whose coordinate is equal to the y-coordinate 
of that horizontal segment of g in R 2 . 

To apply the WSF algorithms to the WkCenter problem, we need a data 
structure for answering queries q(i,j) = cx(i,j) with 1 < i < j < n. Suppose 
such a data structure can be built in 0(n(n)) time and can answer each query 
in 0{q{n)) time; then we say the time bounds of the data structure are 
0(-7r(n), q{n)). The two lemmas below follow from the results in [5]. 

Lemma 2. [3] Suppose there is a data structure for the queries a(i,j) with time 
bounds 0(yr(n), q{n)); then the lDkC enter problem is solvable in 0(jr(n)+n-q(n)) 
time. 

Lemma 3. [5] Suppose there is a data structure for the queries a(i,j) with time 
bounds 0(n(n),q(n)); then the lDkC enter problem is solvable in 0(n(n) +q(n) ■ 
fc 2 log 2 f ) time. 

Refer to [5] for the details of the algorithms in the above two lemmas. By 
the above two lemmas, Lemma 2] follows. 

Lemma 4. Suppose there is a data structure for the queries a(i,j) with time 
bounds 0(n(n), q(n)); then the lDkC enter problem can be solved in 0(min{7r(n) + 
n ■ q(n),ir(n) + q(n) ■ k 2 log 2 ^}) time. 

A data structure based on fractional cascading [4] was given in [5] for an- 
swering the queries a(i,j) with time bounds 0(n log n, log 2 n). Consequently, 
by Lemma |4j the lDkCenter problem is solvable in 0(min{nlog 2 n, nlogn + 
k 2 log 2 j log 2 n}) time. In Section [31 we develop a data structure for processing 
the queries a(i,j) with time bounds 0(n log n, log n), which allows us to solve 
lDkCenter in O(nlogn) time. 

The reason why we can solve the lDkCenter problem faster than simply 
applying the WSF algorithms [5] is that the WSF instance constructed above 
from the problem lDkCenter has a special property: The ^-coordinates of the 
points p'i ,P2, ■ . ■ ,p' n are increasing. As shown in Section [3l this special property 
allows us to design a new data structure for the a(i, j) queries with time bounds 
0(n log n, log n). Note that this special property does not hold for the general 
WSF problem studied in [S]. 

3 The Data Structure for Computing az(i,j) 

In this section, we present a data structure with time bounds 0(n log n, log n) 
for answering the a(i,j) queries. In the following, we first model the problem 
of computing a(i,j) as the problem of finding a lowest point in the common 
intersection of a set of half-planes (i.e., the 2-D sublist LP query). 
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Fig. 2. Illustrating the common intersection of the half-planes defined by three points 
pi, P2, and P3. The point p* is the lowest point in the common intersection. 

3.1 The Problem Modeling 

Consider a point subset Pij C P with i < j. Recall that L(pi) < L(pi + i) < 
■■■ < L(pj). To compute a(i,j), we need to find a point / such that the value 
ip(Pij,f) = maxi<t<j(w(pt) ■ d(f,p t )) = maxi< t < j(w(p t ) ■ \L(f) - L(p t )\) is 
minimized and oe(i,j) is the minimized value. Consider an arbitrary point /' on 
L. Since i/j (P^, f) = m&Xi<t<j(w(pt)-\L(f')-L(p t )\), each point p t G P io defines 
two constraints: w(p t ) ■ (L(f) - L(p t )) < ip(P l3 ,f) and -w(p t )-{L(f')- L{p t )) < 

Consider a 2-D icy-coordinate system with L as the x-axis. For each point 
Pt G P, the inequality w(pt) ■ \x — L(jp t )\ < y defines two (upper) half-planes: 
w(p t ) ■ (x — L(p t )) < y and —w(p t ) ■ (x — L(p t )) < y. Note that the two lines 
bounding the two half-planes intersect at the point p t on L. 

Based on the above discussion, if p* — (x* , y*) is a lowest point in the common 
intersection of the 2(i — j + 1) (upper) half-planes defined by the points in Py, 
then a(i,j) — y* and L(f) = x* is the coordinate of an optimal center / on L 
for P^ . Figure [5] shows an example in which each "cone" is the intersection of 
the two upper half-planes denned by a point in P^. Clearly, this is an instance 
of the 2-D linear programming (LP) problem, which is solvable in 0(j — i + 1) 
time |9ll6j . However, we can make the computation faster by preprocessing. Let 
Hp = {hi, fi2, . . . , ti2n} be the set of 2n (upper) half-planes defined by the n 
points in P, such that for each 1 < i < n, the demand point pi defines ft.2i-i and 
hn- Then to compute a(i,j), it suffices to find the lowest point p* in the common 
intersection of the half-planes defined by the points in P^, i.e., the half-planes 

in i?2i-l,2j = {h>2i-l> h>2i, ^2i+l, ^2i+2) • • • , ^2j-l, ^2j}- 

We actually consider a more general problem: Given in the plane a set of n 
upper half-planes H = {hi, h,2, • ■ • , h n }, each query q(i,j) with 1 < i < j < n 
asks for a lowest point p* in the common intersection of all half-planes in Hij = 
{h t i < t < j} C H. We call this problem the 2-D sublist LP query. Based 
on our discussion above, if we solve the 2-D sublist LP query problem, then the 
a(i, j) queries for the lDkCenter problem can be processed as well in the same 
time bound. A data structure for the 2-D sublist LP query problem with time 
bounds 0(nlogn, log 2 n) was given in [5]. 

Yet, the 2-D sublist LP query problem for lDkCenter is special in the follow- 
ing sense. For each half-plane h G Hp for lDkCenter, we call the x-coordinate 
of the intersection point between L (i.e., the a:-axis) and the line bounding h the 
x-intercept of h (or its bounding line). As discussed above, for each point p t G P, 
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the X- intercepts of both the half-planes h^t-i and hn defined by pt are exactly 
the point pt . Since all points of P are ordered along L from left to right by their 
indices, a special property of Hp is that the x-intercepts of all half-planes in Hp 
are ordered from left to right on L by the indices of the half-planes. For a set H 
of half-planes for a 2-D sublist LP query problem instance, if H has the above 
special property, then we say that H is x-intercept ordered. 

Below, we show that if H is x-intercept ordered, then there is a data structure 
for the specific 2-D sublist LP query problem with time bounds 0(n log n, log n). 
Henceforth, we assume that H is an x-intercept ordered half-plane set. In the 
lDkCenter problem, since all the point weights for P are positive Unite values, 
the bounding line of each half-plane in Hp is neither horizontal nor vertical. 
Thus, we also assume that no bounding line of any half-plane in H is horizontal 
or vertical. Again, this assumption is only for simplicity of discussion. 

In the next section, we solve the 2-D sublist LP queries by reducing it to 
computing the convex hull of a query sub-path of a given simple path [13] . Given 
a simple path in the plane, based on compact interval trees, data structures are 
proposed in [13 to compute (in logarithmic time) the convex hull of a query 
subpath that is specified by the indices of the beginning vertex and the end 
vertex of the subpath, and the convex hull is represented (by a compact interval 
tree) such that standard convex hull queries on it can be done in O(logn) time. 

3.2 Answering 2-D Sublist LP Queries 

For each half-plane hi G H, we denote by I (hi) the bounding line of hi] let l(H') 
be the set of the bounding lines of the half-planes in any subset H 1 C H. 

Our problem reduction utilizes a duality transformation [5] , which is a tech- 
nique commonly used in computational geometry, as follows. Suppose we have a 
primary plane V . For each point (a, 6) £ V, it corresponds to a line y — ax — b 
in the dual plane V; the line is also called the dual of the point and vice versa. 
Similarly, each line y = a'x — b' in V corresponds to a point (a' , b 1 ) in V. 

Suppose all half-planes in H are in the primary plane V . By duality, we can 
obtain a set H* of points in the dual plane V corresponding to the lines in 1(H). 
For each query q(i,j) on H, our goal is to locate the lowest point p* in the 
common intersection of all half-planes in Hij (note that due to our assumption 
that no line in 1(H) is horizontal, there is only one lowest point in the common 
intersection). Since all half-planes in H are upper half-planes, an observation 
is that p* is also the lowest point of the upper envelope of the arrangement of 
the lines in l(Hij). Denote by Uij the above upper envelope for l(Hij). Denote 
by H*j the set of points (in the dual plane V) dual to the lines in i(ify). Let 
C*j denote the lower hull of the convex hull of Hfj. It is commonly known [8] 
that the dual of Uij is exactly C*j (in the dual plane T>). Therefore, if we have 
a representation of C?- that can support standard binary-search-based queries, 
then we can compute the lowest point p* in logarithmic time accordingly. 

One may attempt to design a data structure for querying the lower hull on 
any subset HL of H*. However, there are difficulties when doing so "directly", 
which will be explained later. Instead, we use an "indirect" approach, as follows. 
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Recall that we have assumed the bounding line of any half-plane in H is not 
vertical or horizontal. We partition the half-plane set H into two subsets Hi 
and H2 such that a half-plane hi of H is in Hi (resp., H2) if and only if the 
slope of l(hi) is negative (resp., positive). Accordingly, for each subset Hij, we 
also have Hh and Hfj, and we define the envelopes U~ and Ufj accordingly. 
Since the bounding lines of all half-planes in H^ have negative slopes, the upper 
envelope JZy is monotone decreasing from left to right. Similarly, the upper 
envelope Ufj is monotone increasing from left to right. Hence, it is easy to see 
that the lowest point p* is the single intersection of and Ufj. Let H*l and 
H*? be the sets of points in the dual plane corresponding to the lines in i(-ff^) 
and l{Hh), respectively. Denote by G*j and C*f the lower hulls of H*j and H*^, 
respectively. By duality, the intersection of Uj, and {/?■ corresponds exactly to 
the common tangent line of the two lower hulls C*j and C*^ such that both hulls 
are above the tangent line. Note that since all lines in H^ have negative slopes, 
by duality, all points in H*j- are to the left of the y-axis in the dual plane, and 
thus the lower hull Cy is to the left of the y-axis. Similarly, the lower hull C*^ 
is to the right of the y-axis. Namely, the two lower hulls C** and C*^ are on 
different sides of the y-axis. This property can make our computation of their 
tangent line easier. 

In summary, if we can represent both Cy and C*? in such a way that the 
common tangent line can be found efficiently, then p* can be obtained immedi- 
ately. Our remaining task is to derive a way to support convex hull (or lower 
hull) queries on any subset of consecutive points in iJj* (and similarly on 
Our result is that a data structure can be built in 0(n log n) time such that given 
any i < j, the lower hull C*J can be obtained (implicitly) in O(logn) time and 
is represented in a way that supports binary-search-based queries (e.g., compute 
the common tangent of it and another lower hull, say C*f). The details are given 
below. 

Without loss of generality, we assume the bounding lines of the half-planes 
in H all have negative slopes (i.e., H = H\) and the other case can be handled 
analogously. Suppose the bounding line of each hi £ H corresponds to the point 
h* € H* in the dual plane V. Let 7 be the path by connecting all pairs of 
two consecutive points in H* by line segments, i.e., connecting h* to h* +1 for 
i = 1, 2, . . . , n — 1. Consider the line segment connecting h* and h* +1 and the 
line segment connecting h* and h* +1 ; then the two segments are adjacent to each 
other if i + 1 = j or j + 1 = i. The following Lemma [5] shows that the path 7 is 
a simple path, that is, any two line segments of 7 that are not adjacent do not 
intersect. As can be seen from the proof of Lemma[3 we note that the correctness 
of Lemma [5] heavily relies on two properties of H: (1) H is x-intercept ordered; 
(2) the slopes of the bounding lines of the half-planes in H are all negative (or 
positive). Without either property above, the lemma would not hold, and the 
second property also explains why we need to partition the original set H into 
Hi and H 2 - 

Lemma 5. The path 7 is a simple path. 
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Fig. 3. Illustrating the two double edges dw(si) and dw(sj) (the shaded regions) in 
the primary plane V corresponding to the two segments and Sj in the dual plane T>. 
The two segments Xi and Xj are their intersections with the a:-axis. 

Proof: Consider two segments Sj and Sj where Sj connects h* and h* +1 and Sj 
connects h* and h* +1 . Suppose Si and are not adjacent. To prove the lemma, it 
is sufficient to show Si and Sj does not intersect. Since Si and Sj are not adjacent, 
either i + l<joij + l<i. Without loss of generality, assume i + 1 < j. 

Note that Sj and Sj are in the dual plane T>. It is commonly known [S] that 
the dual of in the primary plane V is the double wedge bounded by the lines 
l(hi) and l(hi-\-i) such that the double wedge does not contain a vertical line 
(e.g., see Fig. [3J ; we denote the double wedge by dw(si). Similarly, the dual of 
Sj is the double wedge dw(sj) bounded by l(hj) and l(hj + i). 

Assume to the contrary the two segments Si and Sj intersect each other, say, 
at a point p. Then, p corresponds to a line in the primary plane V, which is in the 
common intersection of the two double wedges dw(si) and dw(sj). However, we 
claim that the common intersection of dw(-Si) and dw(sj) does not contain any 
line, which incurs contradiction. Hence, Si cannot intersect Sj and the lemma 
follows. Below, we prove the above claim. 

Suppose to the contrary there is a line / contained in dw(si)r\dw(sj). Denote 
by Xi (resp., Xj) the intersection of the x-axis and dw(si) (resp., dw(sj)), e.g., 
see Fig. [3] Hence, both Xi and Xj are line segments on the x-axis. Since the 
intersection of l(hj) and the x-axis is strictly to the right of the intersection of 
l(hi + i) and the x-axis (since H is x- intercept ordered), Xi does not intersect Xj. 
Since the slopes of both I (hi) and l(hi+i) are negative, a line contained in dw(si) 
must intersect x,, and thus the line I intersects x, (due to I C dw(si) D dw(sj)). 
Similarly, I also intersects Xj. Hence, we obtain that / intersects both Xj and Xj. 
Since both Xi and Xj lie in x-axis and Xi does not intersect Xj, the line I has to 
be the x-axis. However, since the slopes of both l(hi) and l(h i+ i) are negative, 
the double wedge dw(si) cannot contain the x-axis and thus cannot contain I. 
Therefore, we obtain contradiction and the claim follows. □ 

In light of Lemma El we can utilize the results in [13]. Given a simple path in 
the plane, compact interval tree data structures are proposed in [13] to compute 
the convex hull of a query subpath that is specified by the indices of the beginning 
vertex and the end vertex of the subpath. If applied to 7 in our problem, then 
after spending 0(n log n) time sorting the points in H* by their x-coordinates, we 
have the following results: A data structure can be constructed in 0(n log log n) 
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time that can compute the lower hull Ch for any query q(i,j) in O(logn) time 
and C* is represented (by a compact interval tree) such that any standard 
convex hull queries on C*j can be done in O(logn) time, where the standard 
convex hull queries includes point-in-polygon tests, finding intersections with 
lines, finding tangents through query points, finding extreme vertices in a query 
direction, detecting intersections of two polygons, and finding common tangents 
of two convex hulls. Further, another data structure of construction time O(n) 
and query time O(lognloglogn) is also given in |13) ; in addition, there is also a 
data structure of construction time 0(n log* n) and query time 0(log n log* n) by 
making trade-off between the construction and query [13^. Both data structures 
are applicable to our problem. 

By duality, the x-coordinate of each point h* £ H* corresponds to the slope 
of the line I (hi). Thus, a sorted order of the points in H* by x-coordinate corre- 
sponds to a sorted order of the bounding lines of the half-planes in H by slope. 
The following lemma summarizes our discussions above. 

Lemma 6. In 0(n\ogn) time, we can build a data structure that can answer 
each 2-D sublist LP query in O(logn) time. Further, if the bounding lines of 
the half-planes in H are sorted by their slopes, then there exist three data struc- 
tures for the 2-D sublist LP queries, whose construction time complexities are 
0(n), 0(n\og* n), and O(nloglogn), respectively, and query time complexities 
are Oilogn log log n), 0(lognlog* n), and 0(log n), respectively. 

According to the reduction procedure from the IDkCenter problem to the 
planar points approximation problem [6] , the slopes of the bounding lines of the 
half-planes in H correspond to the weights of the points in the input point set 
P. Therefore, we have the following corollary. 

Corollary 1. There exists a data structure of time complexity 0(n log n, log n) 

for the a(i,j) queries. Further, if the points ofV are sorted on the line L and the 
weights of the points in V are also sorted, then we have three data structures for 
thea(i,j) queries of time complexities 0(n, log n log log n), 0(n log* n, logn log* n), 
and 0(n log logn, logn), respectively. 

By Lemma 21 we have the following result for the IDkCenter problem. 

Theorem 1. The IDkCenter problem is solvable in O(nlogn) time. Further, if 
the points ofV are sorted on L and the weights of the points in V are also sorted, 
then the IDkCenter is solvable in 0(min{n + k 2 log 2 ? log n log log n, n log* n + 
k 2 log 2 ? log n log* n, n log log n + k 2 log 2 ? log n}) time. 

Therefore, if the points of V are sorted on the line L and the weights of the 
points in V are also sorted, for small k (e.g., k = O ( ( log d - 1 " t - ) 1 / 2 ) , which is 
true in many applications), the IDkCenter problem is solvable in 0(n) time. 
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4 Conclusion 



In this paper, we give an 0(n\ogn) time algorithm for the k-center problem on 
a real line. In certain special cases, we can solve the problem in linear time. 

As suggested by Tamir [19], when k = n — 1, the discrete unweighted Wk- 
Center is equivalent to the Min Gap problem, i.e., finding the closest pair of 
neighbors in P. Hence, there is an O(nlogn) time lower bound on the discrete 
version. In fact, by the reduction from the Min Gap problem, we can also show 
that the continuous unweighted lDkCenter also has an 0(n\ogn) lower bound 
on the running time. Indeed, in any optimal solution OPT for k = n — 1, there 
must be a center at the middle position of the closest pair of neighbors in P, and 
that center serves both neighbors; further, any other center in OPT serves one 
and only one demand point in P. Therefore, given an optimal solution OPT, 
since the demand points served by each center are known, the two demand points 
served by the same center are the closest neighbors in P. We thus obtain the 
0(nlogn) time lower bound on the continuous unweighted lDkCenter. 

Note that the linear time algorithms for the unweighted continuous/discrete 
k-center problem on trees [TT] do not violate the 0{n\ogn) time lower bound 
discussed above because the tree structure already gives a partial order of the 
nodes in the tree. An open problem is whether the techniques given in this paper 
can be extended to the tree structure. 
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